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ABSTRACT 

We use data drawn from the DEEP2 Galaxy Redshift Survey to investigate the relationship be- 
tween local galaxy density, stellar mass, and rest-frame galaxy color. At z '^ 0.9, we find that the 
shape of the stellar mass function at the high-mass (logio(M*//i~^ ^q) > 10.1) end depends on the 
local environment, with high-density regions favoring more massive systems. Accounting for this stel- 
lar mass-environment relation (i.e., working at fixed stellar mass), we find a significant color-density 
relation for galaxies with 10.6 < logio(M*//i~^ M©) < 11.1 and 0.75 < z < 0.95. This result is 
shown to be robust to variations in the sample selection and to extend to even lower masses (down 
to log]^Q(M,//i~^ M©) ~ 10.4). We conclude by discussing our results in comparison to recent works 
in the fiterature, which report no significant correlation between galaxy properties and environment 
at fixed stellar mass for the same redshift and stellar mass domain. The non-detection of environ- 
mental dependence found in other data sets is largely attributable to their smaller samples size and 
lower sampling density, as well as systematic effects such as inaccurate redshifts and biased analysis 
techniques. Ultimately, our results based on DEEP2 data illustrate that the evolutionary state of a 
galaxy at z '^ 1 is not exclusively determined by the stellar mass of the galaxy. Instead, we show that 
local environment appears to play a distinct role in the transformation of galaxy properties at z > 1. 
Subject headings: galaxies:statistics, large-scale structure of universe 



1. INTRODUCTION 

With the emergence of large spectroscopic surveys at 
intermediate r edshift, such as the DEEP2 Galax y Red - 
shift Survey (jDavis et all 120031: iNewman et"all I2OIOD 
the y iMOS VLT Deep Survey (WPS. ILe Fevre'eTaLl 
[2OO5I) . and zCOSMOS (|Lillv et all I2007D . robust svs- 
tematic studies of galaxy clustering are now able to be 
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extended out to z ^ 1 (e.g., ICoil et all I2004U 120061: 
IMeneux et al.|[2008ll2009D . In particular, these unprece- 
dented data sets have enabled the local density of galax- 
ies (i.e., on scales of ^ 1-2 h~^ Mpc) to be statisti- 
cally measured over a broad and continuous range of 
environments extending from voids to rich groups and 
poor clusters at z ^ 1 (e.g.. ICooper et al.l 120051 120061 : 
ICucciati et"alll2006l: IKovac et al.ll2010[) . Moreover, indi- 
vidual galaxy groups at z ~ 1 are now able to be reliably 
identified in redshift space, enabling the "field" popula- 
tion to be studie d in relation to the "group" or "cluster" 
population fe g.. IGerk e et al.l [200l Knobel etHI 120091: 
ICucciati et al.ll2009D . 

Capitalizing on the ability to characterize environment 
over half of the age of the Universe, many studies have 
addressed the role of environment in galaxy evolution 
at z < 1 by studying the relationship between environ- 
ment and galaxy p roperties such as r e st-frame color and 
morphology (e.g JGCTke et al.l 120071: lElbaz etHI 120071: 
ICapak etal1l2007l: ICooper et al.ll20d8l ). Fo r example, us- 
ing da ta drawn from the DEEP2 survey, ICooper et al.l 
()2006i) showed that all features of the global correlation 
between galaxy color and environment measured locall y 
are already in place at z '-^ 1 (see also ICoil et al.|[2008[) . 
This work emphasized that physical processes specific to 
clusters (such as ram-pressure stripping and harassment) 
are not required to explain the color-density relation at 
z ~ 1, given the lack of such extreme environments in 
the DEEP2 sample. 

However, as highlighted by several recent analyses, 
nearly all of the early studies addressing the relation- 
ships between environment and galaxy properties at in- 
termediate redshift primarily utilized luminosity-selected 
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samples and not stellar mass-selected samples. A num- 
ber of recent papers have concluded that while there is 
evidence for a color-density relation at fixed luminosity 
at 2 '^ 1, they fin d no color-densi t y rela tion at fixed 
mass. For example. iScodeggio et al.l (|2009( ) employ data 
from VVDS to examine the relationship between rest- 
frame B — I color and environment as a function of 
redshift within fixed bins of stellar mass. They find no 
evidence for a color-density relation at fixed mass over 
the entire redshift (0.2 < z < 1.4) and stellar mass 
(9 < logig(M*) < fl) ranges probed by their data. 
Using zCOSMOS to characterize the local environment, 
iTasca et all ()2009[) also find that at 0.5 < z < 1 there 
is no variation in galaxy morphology (i.e., early- versus 
late-type) as a function of local galaxy density at stellar 
masses of logjg(M*//i~^ Mq) > 10.5. From these re- 
sults, both of the aforementioned studies conclude that 
the properties (i.e., color and morphology) of massive 
galaxies are independent of environment at 2; ~ 1 . 

Given the substantial noise in all current measures 
of environment (both statistical measures of the local 
galaxy density and group/field catalogs), intrinsic corre- 
lations between galaxy properties and environment can 
easily be smeared out, such that no significant trends 
with environment are detectable in the data. That is, 
it is far easier to smear out a correlation with environ- 
ment than it is to erroneously detect one as, generally, 
measures of local density are dependent on redshift mea- 
surements of neighboring galaxies, while a given galaxy's 
properties are measured from its own imaging and spec- 
troscopic data. It is difficult to imagine a systematic 
effect that would yield a false correlation between inde- 
pendently determined quantities. 

Conclusions based on the lack of an apparent envi- 
ronment dependence must therefore be drawn with cau- 
tion. In this paper, we use data from the DEEP2 survey 
to investigate the significance or robustness of the re- 
cent, aforementioned conclusions inferred from studies of 
galaxy environments at intermediate redshift. In Section 
[21 we describe our data set, with results and discussion 
presented in Sections [3] and 21 respectively. Throughout, 
we employ a ACMD cosmology with w = —1, Q.m = 0.3, 
f^A = 0.7, and a Hubble parameter of Hq = 100 /i km 
s^^ Mpc~^. All magnitudes are on the AB system. 

2. DATA 

Among the current generation of deep spectroscopic 
redshift surveys at z ~ 1, the DEEP2 Galaxy Redshift 
Survey provides the largest sample of accurate spectro- 
scopic redshifts, the highest-precision velocity informa- 
tion, and the highest sampling densitvF^ Altogether, 
these attributes make the DEEP2 survey the best ex- 
isting spectroscopic data set with which to characterize 
local environment. In this paper, we utilize a sample of 
23,037 galaxies wit h accurate redshifts ( quality Q = 2> 
or 4 as defined by iNewman et al.l |2010[ ) in the range 
0.75 < z < 1.25 and drawn from aU four of the DEEP2 
survey fields. 

For each galaxy in the DEEP2 sample, rest-frame 
U ~ B colors and absolute _B-band magnitudes, M^, 

^^ Note that the sampUng density for a survey is defined to be 
the number of galaxies with an accurate redshift measurement per 
unit of comoving volume and not the number of galaxies targeted 
down to an arbitrary magnitude limit. 



are ca lculated from CFHT BRI photometry (jCoil et al.l 
2004al) using t he j^ -correction procedure described in 



Willmer et all (|2006[ ). For a portion of the DEEP2 cat- 



alog, stellar masses may be calculated by fitting spectral 
energy distributions (SEDs) to WIRC/Palomar J- and 
A's-band photometry in conjunction with the DEEP2 
BRI data, a c cording to the prescriptions described by 
iBundv et al.l (120051 . I2006D . However, the near- infrared 
photometry, collected as part of the Palomar Observa - 
tory Wide-field Infrared fPOWIR. IConselice et al.ll200l 
survey, does not cover the entire DEEP2 survey area, 
and often faint blue galaxies at the high-z end of the 
DEEP2 redshift range are not detected in Kg . Because 
of the se two effects, the stellar masses of IBundv et al.l 
()2006[) have been used to calibrate stellar mass estimates 
for the full DEEP2 sample that are based on rest-frame 
Mb and B ~ V values derived from the DEEP2 dat a 
in conjunction with the expressions of iBell et al.l (J2003D . 
We empirically correct these stellar mass estimates to 
the IBundv et al.l ()2006D measurements by accounting fo r 
a mild color and redshift dependence ()Lin et al.ll2007D : 
where they overlap, the two stellar masses have an rms 
difference of approximately 0.3 dex after this calibration. 

To characterize the local environment, we compute 
the projected third-nearest-neighbor surface density (E3) 
about each galaxy in the DEEP2 sample, where the sur- 
face density depends on the projected distance to the 
third-nearest neighbor, £'p,3, as E3 = 3/(7rZ?p 3). In com- 
puting E3, a velocity window of ±1000 km s~^ is utilized 
to exclude foreground and background galaxies along the 
line of sight. Varying the width of this velocity win- 
dow (e.g., using ±1500 km s~^) or tracing environment 
according to the projected distance to the fifth-nearest 
neighb or has no significant effect on our results. In the 
tests of ICooper et al.l (|2005[ ). this projected n*''-nearest- 
neighbor environment estimator proved to be the most 
robust indicator of local galaxy density for the DEEP2 
survey. 

To correct for the redshift dependence of the DEEP2 
sampling rate, each surface density value is divided by 
the median E3 of galaxies at that redshift within a win- 
dow of Az = 0.04; correcting the measured surface den- 
sities in this manner converts the E3 values into mea- 
sures of overdensity relative to the median density (given 
by the notation 1 + <53 here) and effectivel y accounts for 
the re dshift variations in the selection rate (jCooper et al.l 
l2005f l. Finally, to minimize the effects of edges and 
holes in the survey geometry, we exclude all galaxies 
within \ h~^ comoving Mpc of a survey boundary, reduc- 
ing our sample to 17, 767 galaxies in the redshift range 
0.75 < z < 1.25. 

To investigate the relationship between galaxy prop- 
erties and environment within the high-mass segment 
of the galaxy population, we define a subsample 
of galaxies with stellar mass in the range 10.6 < 
logio(M,//i-2 Mq) < 11.1 and a redshift of 0.75 < 
z < 1.05. The median redshift for this subsample of 
1586 galaxies is 0.89 and the median stellar mass is 
logjQ(M,//i~^ M0) ~ 10.8. We restrict the subsample 
in redshift to a range over which the DEEP2 redshift 
selection function is relatively flat. However, over this 
redshift range the sample is incomplete at the adopted 
mass limit. For example, at z = 0.9 the Rab = 24.1 mag- 
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Fig. 1. — the rcst-framc U — B versus Mb color-magnitude dis- 
tribution for DEEP2 galaxies in the spectroscopic sample within 
the redshift range 0.75 < z < 1.25. The three solid black vertical 
linos show the completeness limit of the survey at z = 0.9, 1.05, 
and 1.25, while the solid and dashed red lines show linos of con- 
stant stellar mass corresponding to log]^Q{M*/h~'^ Mq) = 10.6 and 
10.8, respectively. The dashed blue/red horizontal line shows the 
division between the red sequence and the blue cloud as given by 
Equation 19 from lWillmer et al.l 1I2OO6I '). 

nitudc limit of DEEP2 includes all galaxies with stellar 
mass > 10^"'^ M»//i^^ M© independent of color, but pref- 
erentially misses red galaxies at lower masses (of. Figure 

m. 

Figure [1] shows the distribution of galaxies in color- 
magnitude space for the entire DEEP2 sample at 0.75 < 
z < 1.25, with lines of constant stellar mass overlaid. 
As sh own in many previous studies (e.g.. | Bell fc de Jong 
[200lHHogg et al.(l2003l: ICooper et al1l2008[ ). selecting by 
stellar mass as opposed to optical luminosity selects a 
different portion of the galaxy population due to the 
dependence of stellar mass-to-light ratio on rest-frame 
color. A i?-band luminosity-selected sample is biased to- 
wards lower stellar masses for blue galaxies and higher 
stellar masses for red galaxies. 

Now, both locally and at z ~ 1, the local environ- 
ment of galaxies on the red sequence has been shown to 
depend on luminosity such that more lur ninous systems 
favor higher-density regi ons on average ([Blanton et al.l 
120051 : iCooper et al.l 120061 ). Thus, relative to a random 
population of blue galaxies, red galaxies of equivalent 
stellar mass will be (less luminous and) in lower-density 
environments than red galaxies selected to have the same 
luminosity. That is, taking into account the known cor- 
relation between local environment and galaxy luminos- 
ity, for a galaxy sample selected to have fixed stellar 
mass (versus fixed luminosity) one expects to find a 
weaker color-density relation. However, it remains un- 
clear whether or not the correlation between galaxy en- 
vironment aiidj^st^framecolor found at fixed luminosity 
at 2; ~ 1 ([Cooper et al.|[2006() is entirely a projection of 
the correlation between environment and luminosity in 
concert with the dependence of stellar mass-to-light ra- 
tio on color such that at fixed stellar mass there is no 
color-density relation. As discussed in EJH several recent 
analyses have arrived at this conclusion. 

3. ANALYSIS 

As highlighted in fjl] the goal of this work is to study 
the relationship between galaxy color and local environ- 
ment at fixed stellar mass. Recent analyses using a vari- 



ety of data sets have approached this task using galaxy 
samples in stellar-mass-selected bins across a range of 
environments. Here, we first investigate whether this 
methodology is appropriate. Finding it is not — that 
is, finding that the stellar mass function is different in 
different environments — we then apply improved tech- 
niques which account for the correlation between stellar 
mass and environment in an examination of the color- 
density relation at z -^ 1. In contrast to recent results 
in the literature, we find that there is a significant color- 
density relation at fixed stellar mass. The most likely 
reason for this discrepancy is that both larger measure- 
ment errors and systematics in other data sets smear 
out the underlying environment dependence. Finally, we 
emphasize that absence of evidence is not evidence of ab- 
sencejill that is, not finding a correlation between color 
and environment at fixed stellar mass does not directly 
indicate an intrinsic lack of correlation. By making such 
an assumption, the conclusions of several recent analyses 
have been rendered difficult to interpret. 

3.1. The Relationship between Stellar Mass and 
Environment 

In order to study the relationship between galaxy prop- 
erties and environment at fixed stellar mass, the galaxy 
sample under study is often restricted to a narrow range 
in stellar mass such that correlations between stellar 
mass and environment are negligible. However, at in- 
termediate redshift, sample sizes are generally limited 
in number such that using a particularly narrow stellar 
mass range (e.g., ~0.1-0.2 dex in width) significantly re- 
duces the statistical power of the sample. For this reason, 
broader stellar mass bins (e.g., ^-^ 0-5 dex) are commonly 
employed (e.g..lScodeggio et al.ll2009Hlovino et al.ll2009i: 
iPatel et al.ll2009HKovac et al.ll2009l : iMaltbv et al.ll2010( ). 
Now, if the shape of the stellar mass function depends 
on environment, then the typical stellar mass within a 
broad mass bin may differ significantly from one density 
regime to another. Such an effect would clearly impact 
the ability to study the color-density relation at fixed 
stellar mass. 

To investigate the level to which we are sensitive to 
a stellar mass-environment relation within our adopted 
stellar mass bin, we select those galaxies within the top 
10% of the overdensity distribution for all galaxies at 
10.6 < logio(M*//i-2 Mq) < 11.1 and 0.75 < z < 1.05 
— a subsample of 159 galaxies. From the correspond- 
ing bottom 50% of the overdensity distribution, we draw 
1000 random galaxy subsamples (each consisting of 159 
galaxies) , where each subsample is selected to match the 
redshift distribution of the galaxies in the high-density 
subsample. 

By matching in redshift, we remove the projection of 
any possible residual correlation between our environ- 
ment measurements and redshift (in concert with the 
known redshift dependence of the survey's stellar-mass 
limit) onto the observed stellar mass- environment rela- 
tion. This effect ap pears to be impacting the analysis of 
llovino et al.l ()2009[ ). As shown in Figure 12 of that work, 

^^ We note that this turn of phrase, which we utilize in the 
title to this worlt, is a commonplace made popular by Carl Sagan 
among others in discussing a particular logical fallacy related to 
the misinterpretation of non-significant scientific results. 
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the typical redshift of "isolated" galaxies (i.e., galaxies in 
low-density environments) is systematically skewed to- 
wards higher z relative to the sample of "group" mem- 
bers (i.e., galaxies in high-density environments). This 
effect is likely due to the decreasing sampling density of 
a magnitude- limited redshift survey at higher z, which 
introduces a bias against identifying groups at higher 
redshift. Note that our adopted stellar m ass bins arc 
well- matched to those of llovino et al.l ()2009D , taking into 
account differences in the adopted value of the Hubble 
parameter. 

To test whether our high-density subsample and the 
random low-density subsamples are drawn from the 
same underlying stellar mass distribution, we apply 
the one-sided Wilcoxon- Mann- Whitney (WMW) U test 
(|Mann fc WhitnevI 119471 ). a non-parametric test that is 
highly robust to non-Gaussianity because it relies on 
ranks rather than observed values. The result of the 
WMW U test is the probability (Pu) that a value of 
the U statistic equal to the observed value or more ex- 
treme would result if the "null" hypothesis (that both 
samples are drawn from identical distributions) holds. 
The WMW U test is particularly useful for small data 
sets (e.g., compared to ot her related tests such as the 
chi-square two-sample test, IWall fc Jenkinj|2003[ ). as we 
have when selecting galaxies from a narrow stellar mass 
range and in extreme environments, due to its insensitiv- 
ity to outlying data points, its avoidance of binning, and 
its high efficiency. Note that since this test is one-sided, 
possible Pjj values range from to 0.5; for a Pu value 
below 0.025 (corresponding closely to 2a for a Gaussian), 
we can reject the null hypothesis (that the two samples 
have the same distribution) at greater than 95% signifi- 
cance. 

In Figure [21 we plot the cumulative distribution of stel- 
lar masses for the 159 sources in the high-density subsam- 
ple along side that for the 1000 random subsamples (each 
consisting of 159 galaxies) matched in redshift but resid- 
ing in low-density environments. Performing a one-sided 
WMW U test on the stellar mass measurements for the 
low- and high-density populations, we find that the stel- 
lar mass distribution for the galaxies in high-density en- 
vironments is skewed to slightly higher stellar mass, with 
a probability of Pjj ~ 0.09. Meanwhile, the cumulative 
redshift distributions for the low- and high-density sub- 
samples, shown in the inset of Figure [21 are well- matched 
with the WMW U test yielding a Pu > 0.48, thus con- 
firming that the redshift distributions for the two sam- 
ples are indistinguishable and therefore not significantly 
contributing to the difference in the stellar mass distri- 
butions. 

The results of the WMW U test are supported by a 
comparison of the mean stellar masses for the low- and 
high-density subsamples, which are distinct at a > 1.5ct 
level, thereby supporting the conclusion that there is a 
non-negligible stellar mass-environment relation within 
the given mass bin. In addition to directly comparing 
the arithmetic means of the stellar mass distributions, 
we also utilize the Hodges-Lehmann (H-L) estimator of 
the mean, which is given by the median value of the 
mean stellar mass computed over all pa irs of galaxies 
in the sample (jHodges fc Lehma"nnlll963[ ). Like taking 
the median of a distribution, the H-L estimator of the 
mean is robust to outliers, but, unlike the median, yields 



results with scatter (in the Gaussian case) comparable 
to the arithmetic mean. Thus, by using the H-L esti- 
mator of the mean, we gain robustness as in the case 
of the median, but unlike the median, our measurement 
errors are increased by only a few percent. In Figure 
[31 we show the distribution of differences between the 
Hodges-Lchmann estimator of the mean stellar mass for 
the high-density subsample relative to that for the 1000 
low-density subsamples, where the median difference in 
stellar mass is AlogjQ(M,//i~^ M©) ~ 0.019 as illus- 
trated by the dotted vertical line. Within the stellar 
mass range of 10.6 < logiQ(M*//i~^ M©) < 11.1, we find 
(by all of these methods) ^ 1.5cr evidence for a weak stel- 
lar mass-environment relation at z ~ 0.9 such that more 
massive galaxies favor higher-density regions. 

This weak correlation between stellar mass and envi- 
ronment is robust to variations in the stellar mass and 
redshift ranges used to select the galaxy population. For 
example, pushing to slightly lower stellar masses (10.2 < 
logio(M*//i~2 Mq) < 10.7), while stih limiting the sam- 
ple to a range of only 0.5 dex, we again find evidence 
for a weak correlation between between stellar mass and 
environment (see Table [1} . Expanding the range of stel- 
lar masses probed to 10.1 < logio(M*//i~^ M©) < 11.1, 
the stellar mass-density relation is increasingly evident, 
such that we find a highly-significant correlation be- 
tween stellar mass and environment at z ^ 0.9, with 
Pu < 0.01. This result is in general agreement with ex- 
isting studies of the projected galaxy corre lation function 
as a function of stellar mass at z ^ 1 (e.g.. lMeneux et aD 
l2008t iFoucaud et al.l |2010[) , which find weak evidence 
that more massive galaxies are more strongly clustered 
on small scales. Likewise, measurements of the stellar 
mass function at slightly lower redshift {z ~ 0.5) also 
suggest a variation i n the shape of the mass f unction with 
local environment ([Bplzonella et al.l 120091 : iKovac et al.l 
[200llPenget al.ll2010() . 

Thus, when looking to quantify environmental de- 
pendencies at fixed stellar mass, th e use of broad 
stellar mass bins (e.g., those used bv iScodeggio et al 



20091. Ilovin q et al.l 120091 iTasca et aD 120091 IKovac et al 



20091 and iGriitzbauch et al.l I2010D or simple mass- 



limited sample selec t ions fe.g.. Tvari der Wei et al.|[2007l : 
iPannella et al.l l2009t iRettura et al.ll2010[ ) is inappropri- 
ate as both approaches implicitly assume that there is 
no variation in the stellar mass function with environ- 
ment. In the following subsection, we develop improved 
techniques that account for the apparent environmen- 
tal dependence of the shape of the stellar mass func- 
tion, thereby enabling an unbiased analysis of the color- 
density relation at fixed stellar mass. 

3.2. The Color-Density Relation at Fixed Stellar Mass 

The ultimate goal of this work is not to study the 
correlation between stellar mass and environment, but 
rather to investigate the relationship between galaxy 
color and environment at fixed stellar mass. But as 
shown in the above analysis, we must be careful to 
account for the weak correlation between stellar mass 
and overdensity even with stellar mass bins as nar- 
row as 0.5 dex. Thus, we now select those galaxies 
within the top 10% of the overdensity distribution for 
ah galaxies at 10.6 < logio(M*//i~2 j^j^) < m a^d 
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Fig. 2. — the cumulative stellar mass and redshift (see inset) distributions for the 159 DEEP2 galaxies comprising the top 10% of 
the environment distribution within the stellar mass and redshift ranges of 10.6 < logj^Q(M,/h~^ Mq) < 11.1 and 0.75 < z < 1.05 in 
comparison to the corresponding cumulative distributions for the 1000 random galaxy subsamples drawn from the lowest 50% of the same 
environment distribution. As discussed in the text, the low-density subsamples, which are each composed of 159 galaxies, are selected to 
have the same redshift distribution as the high-density population. However, the stellar mass distribution is found to be biased in the 
different environment regimes. 
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Fig. 3. — the distribution of differences between the Hodges-Lehmann (H-L) estimator of the mean stellar mass {left) and mean redshift 
(right) for the high-density subsample relative to the corresponding H-L estimator of the mean for each of the 1000 low-density subsamples. 
The median differences in stellar mass and redshift are denoted by the dotted vertical lines, while the dashed vertical lines and corresponding 
error bars show the difference in the arithmetic means (and uncertainty in that difference) for the low-density and high-density populations. 
We find a small offset in stellar mass of Alog]^Q(M«/h~^ Mq) ~ 0.019, while the difference in mean redshift of the two samples is 
consistent with zero (by construction). Within the broad stellar mass range of 10.6 < logjQ(M»//i~-^ M©) < H-l, we detect a weak stellar 
mass-environment relation at 2 ~ 0.9 such that more massive galaxies favor higher-density regions. 



0.75 < z < 1.05 (the same high-density subsample of 
159 galaxies), and from the corresponding bottom 50% 
of the overdensity distribution we randomly draw 1000 
subsamples (each composed of 159 galaxies) so as to 
match the joint redshift and stellar mass distributions 
of the galaxies in the high-density subsample. Mem- 
bers of the low-density subsample are drawn randomly 
from within a two-dimensional window with dimensions 
of Az2 < 4 ■ 10-4 and Alogio(M,)2 < 5 ■ IQ-^ of a 
randomly-selected object in the high-density subsample. 
Varying the size of this window by factors of a few in 
each dimension has no significant effect on our results. 
Given the random nature of the matching, some objects 
are repeated in the low-density subsamples. However, 
for each subsample of 159 galaxies, '--^90% of the galax- 
ies are unique; requiring all members of a subsam ple to 
be unique would skew the statistics (|Efron|[l981[ ). By 
matching our high- and low-density subsamples in stel- 
lar mass as well as redshift, we are able to effectively 
study the correlation between galaxy properties such as 
color and environment at fixed stellar mass. 

As shown in Figures 2] and [51 we find a significant rela- 
tionship at z ~ 0.9 between rest-frame U ~ B color and 
local galaxy overdensity at fixed stellar mass for galaxies 
with 10.6 < logio(M*//i-2 Mq) < 11.1. The cumulative 
color distributions for the two subsamples drawn from 
separate environment regimes are significantly distinct, 
with the galaxies in high-density environments skewed 
towards redder rest-frame colors. The WMW U test con- 
firms what is apparent in Figure SI yielding Pu < 0.002 
when comparing the color distributions of the high- and 
low-density subsamples: these color distributions are dis- 
tinct at > 3(7 confidence. By construction, the corre- 
sponding stellar mass and redshift distributions are in- 
distinguishable (see inset plots in Fig.|4]), with Pjj ~ 0.47 
and Pjj ^ 0.49, respectively. 



In Figure [SJ we show the distribution of differences 
between the Hodges-Lehmann estimator of the mean U~ 
B color, stellar mass, and redshift for the high-density 
subsample relative to the corresponding estimator of the 
mean for the 1000 low-density subsamples. The median 
differences in stellar mass and redshift are consistent with 
zero, Az- -M0-4±0.009 and Alogio(MJ/i-2 Mq) ~ 
0.003 ± 0.012, while the median offset in rest-frame color 
is A([/ ~B) r^ 0.033 ± 0.013 such that galaxies in high- 
density environs are typically redder in color at fixed 
stellar mass and redshift. Similarly, the arithmetic means 
of the two color distributions arc accordingly found to be 
distinct at >3.5ct (see Tabled]). Altogether, the DEEP2 
data show a robust color-density relation at fixed stellar 
mass at z ^ 0.9 for galaxies in the stellar mass regime of 
10.6 < logio(M,//i-2 Mq) < 11.1. 

To test the robustness of our results to the particular- 
ities of the sample selection, we repeat the analysis de- 
scribed above for several samples spanning varying red- 
shift and stellar mass regimes. For example, broaden- 
ing the redshift range over which we select galaxies to 
0.75 < z < 1.25, thereby increasing the size of the sam- 
ple, we again find a statistically significant relationship 
between rest-frame color and environment within our 
adopted stellar mass bin of 10.6 < logio(M*/ft,-^ Mq) < 
11.1. For the 211 galaxies in the high-density regime 
(again the highest 10% of the overdensity distribution) 
at 0.75 < z < 1.25, the cumulative distribution of 
U — B color is skewed towards redder colors relative to 
the comparison set of galaxies in low-density environ- 
ments, yielding Pu < 0.01 and with the means of the 
two color distributions distinct at a ~ 3(t level. In Ta- 
ble m we list the results from similar analyses of sev- 
eral other galaxy samples. When varying the redshift 
and/or stellar mass regimes probed, we continue to find 
a significant color-density relation at fixed stellar mass 
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Fig. 4. — the cumulative rest-frame U — B color, stellar mass, and redshift distributions for the 159 DEEP2 galaxies with the highest mea- 
sured overdcnsitics (the top 10% of the environment distribution) within the stellar mass and redshift ranges of 10.6 < log]^Q{M»/h~^ I^o) < 
11.1 and 0.75 < z < 1.05 in comparison to the corresponding cumulative distributions for the 1000 random galaxy subsamples drawn from 
the lowest 50% of the same environment distribution. As discussed in the text, the low-density subsamples, which are each composed of 159 
galaxies, are selected to have the same stellar mass and redshift distributions as the high-density population (see inset plots). For galaxies 
with stellar masses of 10.6 < log]^Q{M»/h~^ ^q) < H-l at 2 ~ 0.9, we find a significant color-density relation at fixed stellar mass such 
that red galaxies preferentially reside in overdense regions relative to their bluer counterparts. 
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Fig. 5. — the distribution of differences between the Hodgcs-Lehmann (H-L) estimator of the mean U — B color, redshift, and stellar mass 
for the high-density subsample relative to the corresponding H-L estimator of the mean for each of the 1000 low-density subsamples. The 
median differences in U — B, z, and stellar mass are denoted by the dotted vortical lines, while the dashed vertical linos and corresponding 
error bars show the difference in the arithmetic means (and uncertainty in that difference) for the low-density and high-density populations 
(taken as a whole). Recall that the redshift and stellar mass distributions of the low- and high-density samples match by construction. We 
find a significant offset in rest-frame color of A(f/ — _B) ~ 0.04, while the difference in mean redshift and stellar mass for the two samples 
is consistent with zero. 



at 2: ~ 1. This result is found to hold even at slightly 
lower masses; restricting to only those galaxies with 



10.2 < logio(M,//i"2 jvig) < 10.7 (a sample with a me- 
dian stellar mass of log]^Q(M*//i^^ Mq) ~ 10.4), we again 
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find that redder galaxies tend to favor overdense regions 
at fixed stellar mass and redshift. 

4. DISCUSSION 

As shown in Figure IH we find that there exists a cor- 
relation between rest-frame U — B color and local galaxy 
overdensity within the high-mass (logj^Q(M*//i^^ Mq) ~ 
10.8) segment of the galaxy population at z -^ 0.9 such 
that red galaxies favor higher-density environments at 
fixed stellar mass. As noted in Section [TJ several recent 
studies utilizing data from VVDS and zCOSMOS have 
concluded that no such correlation exists within this stel- 
lar mass and redsh ift regime fe.g.. lScodeggio et al.ll2009l : 
llovino et al.l 120091 ) . Here, we discuss the likely reasons 
for the discrepancy between these conclusions and our 
results. 

It is highly unlikely that the highly significant color- 
density relation at fixed stellar mass within the DEEP2 
data is spurious. Rest-frame galaxy colors and stellar 
masses, which depend on the adopted stellar population 
models, photometry, and only coarsely on redshift, are 
entirely independent of the environment measures, which 
depend upon angular position and high-precision redshift 
information (both for the galaxy in question and neigh- 
boring galaxies). There is no reasonable mechanisms 
that would produce a false correlation between these in- 
dependent quantities. For instance, one might suggest 
that contamination of photometry is an issue in dense re- 
gions; however, even in the most overdense environments 
(e.g., the top 5% of the environment distribution), the 
typical distance to the 3'''*-nearest neighbor corresponds 
to ^ 35" on the sky, which is much larger than the aper- 
ture sizes used in photometry. Instead, as highlighted 
in Section [1] it is far more likely that the color-density 
relation apparent in DEEP2 has been smeared out in 
studies using other data sets due to the smaller sample 
sizes employed and the significantly larger errors in the 
environment measures derived from those data sets. 

Many of the limitations regarding the VVDS data 
set, as they relate to detecting correlations between 
galaxy properties and environment at z ~ 1, are dis- 
cussed at length by iCooper et al.l (l2007f) in a com- 
parison to t he work by Cucciati et al.l (I2006D . Using 
DEEP2 data, iCooper et al.l ()2007D fiiida significant cor- 
relation between rest-frame galaxy color and environ- 
ment at Mb < -20.5 - 5 ■ logio(/i) and 0.9 < z < 1.2, 
where non e is evident in the ana lysis of VVDS data pre- 
sented by iCucciati et al.l (|2006( ) . While both of these 
studies employed luminosity-selected galaxy samples, the 
lessons learned from the comparison are no less appli- 
cable to analyses of samples selected according to stel- 
lar mass. He r e, we review the arguments presented by 
ICooper et al.l ()2007l ) as they relate to the more recent 
analyses using VVDS, zCOSMOS, and other data sets. 

At z > 0.75, the sample size of both VVDS and 
zCOSMOS arc significantly smaller than that collected 
by DEEP2. At these redshifts, DEEP2 has more than 
26, 891 secure (95% or greater confidence) redshifts, while 
the VVDS and zCOSMOS data sets include only 1,581 
and 967 high-quality (i.e., fiag:=3,4) redshifts. respec- 
tively. 

In addition to their smaller statistical power relative to 
DEEP2, both VVDS and zCOSMOS have a significantly 
lower number density of tracer objects (cf. , Newman et al.l 



|2010[ ) , which directly influences the level to which the re- 
sulting data sets can trace local galaxy density. Crudely, 
errors on overdensity measurements on some length scale 
should be dominated by Poisson statistics and hence 
should be proportional to the square root of the num- 
ber density of tracers. To test the impact of a re- 
duction in the sampling density, we dilute the DEEP2 
data set by a factor of 3, thereby creating a sample 
that mimics the sampling density of VVDS and zCOS- 
MOS at z ^ 1. Using this tracer population, we then 
recompute the local overdensity about each galaxy in 
the full DEEP2 sample. Note that this approach al- 
lows us to test the role of sampling density indepen- 
dent of sample size. Performing the same analysis de- 
scribed in i )3.2i but using environment measures com- 
puted with the diluted tracer population, we still find 
a color-density relation at fixed stellar mass. However, 
the significance of the trend is noticeably reduced. For 
galaxies with 10.6 < \og^f^{M^ / h'^ Mq) < 11.1 and 
0.75 < z < 1.05, wc find that the color distributions 
for galaxies in high- and low-density environments are 
distinct with a Pjj '^ 0.03 when using the diluted tracer 
population versus Pu < 0.01 when measuring environ- 
ment with the full DEEP2 sample. Comparing the mean 
rest-frame colors as a function of environment, we find 
that the difference in the arithmetic means, A([/ — B), 
between high- and low-density regions is significant at 
only ^ 2cr when employing the diluted tracer sample 
compared to ~ 4a when the full data set is used to define 
environment. Thus, while variation in the sampling den- 
sity may not entirely explain the lack of evidence for a 
correlation between environment and rest-frame galaxy 
color at fixed stellar mass in the VVDS and zCOSMOS 
data sets, it clearly plays a substantial role. As the num- 
ber density of tracers decreases, the uncertainty on each 
environment measure increases such that underlying cor- 
relations with environment are increasingly smeared out. 
Analyses using the VVDS and zCOSMOS data sets 
also often rely on the use of less accurate redshifts. For 
example, the VVDS redshift catalog is domina ted by 
lower-quality redshifts (fiag = 2. lllbert et al.ll2005l ) at z ^ 

1. A number of tests have found that the fiag = 2 red- 
shifts have a non-negli gible error rate (~ 20 %, C. Wolf, 
private communication uLe Fevre et al.ll2005f). The inclu- 
sion o f these sources in the analysis of iScodeggio et all 
(|2009f l results in both inaccurate rest-frame colors (via 
incorrect iiT-corrections) as well as inaccurate overdensity 
measurements in their high-redshift (0.7 < z < 1.4) sam- 
ple and thus contributes to the lack of color-density re- 
lation observed within the 10.2 < login (M^//t~^ Mp^) < 
10.7 stellar mass regime (see Fig. 5 of 'Scodeggio et al.l 

|2009( ). Along the same lines. iTasciTet al. (2009) utilize 
zCOSMOS redshifts with quality codes fiag = 4, 3, 9, and 

2. This includes objects for which only a single emission 
line was identified in the observed spectrum. While these 
spectroscopic redshifts are checked with photometric red- 
shifts and comprise only ~6% of the entire sample, these 
objects play a significantly larger role at z > 0.8, com- 
prising nearly one quarter of the sample at high redshift. 

Given that there are multiple ways to explain the 
smearing out of a real correlation between galaxy color 
and environment, but no reasonable mechanism by which 
to generate a false correlation, we conclude that the 
color-density relation at fixed stellar mass that we ob- 
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TABLE 1 

Summary of Results for Stellar Mass-Density Relation at Fixed Redshift 



Sample 


JVhi 


gh — density 


Puiz) 


Pc(logio{M.)) 


Az 


Alogiu(M.) 




0.75 <z< 1.05 
10.6 < logio(M,) < 11.1 




159 


0.49 


0.09 


0.000 ±0.007 


0.019 ±0.012 


0.75 < z < 1.05 
10.1 < logio(M.) < 11.1 




456 


0.48 


<0.01 


0.000 ± 0.004 


0.051 ±0.012 


0.75 < z< 1.25 
10.6 < logio(M.) < 11.1 




211 


0.49 


0.13 


0.000 ±0.005 


0.012 ±0.010 


0.75 < z < 1.05 
10.2 < logio(M,) < 10.7 




273 


0.49 


0.08 


0.000 ± 0.006 


0.012 ±0.008 





Note. — For a variety of galaxy samples, we list the results of the WMW U test (Pu) and the 
difference in the arithmetic means computed from a comparison of the redshift and stellar mass 
distributions of the respective low- and high-density samples. The P- value, Pu, indicates the prob- 
ability that differences in the distribution of the stated quantity as large as those observed (or 
larger) would occur by chance if the two samples shared identical distributions. The number of 
galaxies in the high-density sample (picked to be the top 10% of the environment distribution) is 
given by A'^high- density Note that the difference in the mean for the galaxy property x is given by 
< 2;high-dcnsity > — < ^^low-density > ^i^d that stellar masses are in units of h~^ Mq. 



TABLE 2 
Summary of Results for Color-Density Relation at Fixed Stellar Mass and Redshift 



Sample 


-^^high — density 


Pu(z) 


Pc/(logio(M.)) 


Pu{U-B) 


Az 


Alogio(M.) 


A(U-B) 




0.75 < z < 1.05 
10.6 < logio(M.) < 11.1 


159 


0.49 


0.47 


<0.01 


0.000 ± 0.007 


0.002 ±0.012 


0.040 ±0.011 


0.75 < z< 0.95 
10.6 < logio(M.) < 11.1 


101 


0.48 


0.47 


0.01 


0.000 ± 0.005 


0.002 ±0.015 


0.043 ±0.014 


0.75 < z< 1.25 
10.6 < logi(,(M.) < 11.1 


211 


0.47 


0.46 


<0.01 


0.000 ± 0.008 


0.002 ±0.010 


0.032 ±0.010 


0.75 < z < 1.05 
10.1 <logi(,(M.) < 11.1 


456 


0.48 


0.46 


<0.01 


0.000 ± 0.004 


0.003 ±0.012 


0.029 ±0.009 


0.85 < 2 < 1.15 
10.1 < logio(M.) < 11.1 


404 


0.47 


0.49 


0.02 


0.000 ± 0.004 


0.002 ±0.013 


0.019 ±0.010 


0.75 <z< 1.05 
10.2 < logio(M.) < 10.7 


273 


0.49 


0.47 


<0.01 


0.000 ± 0.006 


0.000 ± 0.008 


0.032 ±0.011 





Note. — For a variety of galaxy samples, we list the results of the WMW U test (Pu) and the difference in the arithmetic means 
computed from a comparison of the redshift, stellar mass, and color distributions of the respective low- and high-density samples. 
The P- value, Pu, indicates the probability that differences in the distribution of the stated quantity as large as those observed (or 
larger) would occur by chance if the two samples shared identical distributions. The number of galaxies in the high-density sample 
(picked to be the top 10% of the environment distribution) is given by Aiiigii_ density Note that the difference in the mean for the 
galaxy property x is given by < s^high- density > — < 2;io„_density > a-nd that stellar masses are in units of h~'^ Mq. Finally, note 
that the galaxy samples listed here are distinct from those detailed in Table [T] These samples are matched in redshift as well as 
stellar mass. 



serve is genuine. The authors of several recent papers 
appear to have gone too far in concluding that because 
they detect no relation, none exists; this has led to a 
problematic interpretation of their results. 

5. SUMMARY 

Herein, we use data from the DEEP2 Galaxy Red- 
shift Survey to complete a detailed study of the rela- 
tionship between rest-frame galaxy color and local en- 
vironment at fixed stellar mass at intermediate redshift. 
Our principal result is that at fixed stellar mass and red- 
shift we find a significant relationship between rest-frame 
U — B color and local galaxy density at 2; ~ 0.9 within 
the massive (10.6 < logio(M*//i"2 Mq) < 11.1) galaxy 
population. This color-density relation is such that red 
galaxies at a given stellar mass are preferentially found 
in overdense regions, thereby showing that the general 
pattern of environmental dependence of galaxy prop- 
erties at fixed stellar mass in the local Universe (e.g., 
iKauffmann eTal.. ^04,; Baldrv et all [2006t Ivan der Well 



|2008[ ) persists to z '-^ 1. Our findings disagree with some 
recent results, which have found no correlation between 
galaxy properties and environment at fixed stellar mass 
at z > 0.8 (e.g.. iTascaet al.l l2009t llovino et all [20091: 
iScodeggio et alll2009[ ). 

As shown in this work, being unable to detect any 
correlation between galaxy properties and environment 
within a specific data set does not alone indicate a true 
lack of environmental dependence. Uncertainties and 
limitations associated with observations are generally 
such that they smear out the underlying correlation be- 
tween environment and galaxy properties; as a result, 
the observed strength of relationships such as the color- 
density relation are lower limits to the inherent environ- 
mental dependence. In order to glean anything mean- 
ingful from an observed lack of correlation between en- 
vironment and galaxy properties such as color and mor- 
phology, additional tests must be undertaken to illustrate 
that the lack of correlation is inconsistent with being due 
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to observational u ncertainties (e.g.. iCooper et al] 120071 : 
iGerke et al.ll2007t) . Moreover, even an observed decline 
or weakening in the strength of an environmental de- 
pendency with redshift must be treated carefully; the 
observed amplitude of environment trends will weaken 
simply because errors generally increase and sample sizes 
generally decrease with increasing redshift. Much of 
the recent e nvironment- or clustering-related work at 
z ~ 1 ( e.g iCucciati et al.|[20ql: iLe Fevre et al.l [20071: 



[2001" 



lovino et al.' '2009': 'Scodcggi o et al.l 
20091; [Bplzonclla ct al. 2009) omits 
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2009; Ko TOc et al.l 

such analysis, making the associated conclusions difficult 

to interpret. 

Ultimately, in this work, we have shown that the evolu- 
tionary state of a galaxy at z ~ 1 as traced by rest-frame 
U—B color (i.e., specific star-formation rate) is not solely 
determined by the stellar mass of the galaxy. That is, our 
results indicate that environment plays a role (or at the 
least is correlated with a causal physical parameter such 
as dark matter halo mass) in defining the evolutionary 
history of a galaxy at z > 1. As previously stated, ex- 
isting analyses in the local Universe arrive at a similar 
conclusion with regard to the separable roles of stellar 
mass and environment at z > 0. We now know that 
the typically redder colors of galaxies in high-density re- 
gions at z < 1 is not solely attributable to the tendency 
of these extreme environs to host systematically more 
massive galaxies. Instead, the variation in color with 
environment appears to also be partially the result of 



density-dependent effects. Establishing the relative im- 
portance of these two evolutionary drivers (stellar mass 
and environment) is a difficult task, requiring a careful 
examination of the associated measurement errors. Such 
analysis will no doubt be the focus of much future work. 
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